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Abstract 

A susceptibility function is introduced to quantify some aspects of the in- 

termittent stick-slip dynamics of a rough metallic cylinder of length L on a rough 
metallic incline submitted to small controlled perturbations and maintained below 
the angle of repose. This problem is studied from the experimental point of view and 
the observed power-law behavior of x{L) is justified through the use of a general 
class of scaling hypotheses. 
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1 Introduction 



Several dynamic processes of physical [1] , geologic [2] , and technological [3] in- 
terest involve superposed interacting rough surfaces able to display stick-slip 
motion [4]. It has been discovered in the last few years that the intermittent 
sliding or stick-slip dynamics of a rough solid cylinder on a rough inclined sur- 
face submitted to small controlled perturbations is a fluctuation phenomenon 
characterized by non-trivial spatiotemporal scaling laws [5,6], and complex 
critical exponents [7] if the inclination is below the angle of repose. In the 
usual stick-slip dynamics a solid is pulled at a constant driving velocity; in 
the experiments discussed in this paper, on the other hand, the stick-slip dy- 
namics appears as a consequence of a different mechanism: to start the slip 
we resort to small mechanical perturbations on an inclined surface below the 
angle of repose. A solid body on a perturbed incline is an example of a nonequi- 
librium system receiving an incoming energy flow. If energy is continuously 
injected into nonequilibrium systems, a complex sequential response charac- 
terized by time series of events of all sizes is often observed. Besides sliding 
blocks on inclines, other important examples of systems and phenomena asso- 
ciated with a similar type of temporal response include piles of sand and other 
granular materials [8,9]; acoustic emission from volcanic rocks and microfrac- 
turing processes in general [10,11]; interface depinning in magnetic systems 
[12]; stick-slip motion in lubricated systems [13], and turbulence [14], among 
others. 

In this paper, the sliding susceptibility x{L) of a rough metallic cylinder of 
length L placed on a rough inclined surface submitted to small controlled 
perturbations is introduced and examined. To obtain this response function, 
we performed a class of very time consuming experiments in which we count 
the number of perturbations on the incline that give origin to a certain number 
of sliding events of the cylinder. The total number of perturbations in our 
experiments exceeded 56, 000, corresponding to 6, 000 sliding events. Before 
and after each perturbation, the position of the cylinder along the incline 
as well as other variables of interest must be controlled. We find that x{L) 
presents scaling symmetry [15] in a number of physical situations and the 
origin of this behavior is investigated. 

The outline of the article is as follows. The experimental details are described 

in Sec. 2; and in Sec. 3 we report the experimental data, and discuss our 
principal results. Also in Sec. 3, we present scaling hypotheses to explain our 
findings. A summary of our major conclusions and future prospects are found 
in Sec. 4. 
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2 Experimental details 



The basic apparatus used to obtain the experimental data is outhned in Fig. 
1: it consists of a rigid V-shaped anodized aluminum chute made of corner 
plate of 5 mm of thickness and with 90° angular aperture symmetrically dis- 
posed with respect to the vertical plane. The chute was rigidly maintained 
with an inclination 9 with respect to the horizontal, and it is supplied with 
an articulated hammer of mass m — 175g or 75g which hits the base of the 
chute with a controlled (fixed) velocity close to 50cm/s at the moment of the 
impact. The system is mounted on a 200kg table isolated from mechanical 
vibrations. On the chute (with an effective length Lq = 1,500mm) is placed 
a metallic cylinder of length L (5mm to 1,000mm), and the system chute + 
cylinder operates below the critical angle of repose Oc—iaxT^^si where //g is 
the coefficient of static friction of the cyhnder on the chute. Two types of 
chute, hammer, and metallic cylinders were used as will be specified below. 
In the regime of inclination used in the experiments, the series of induced 
sliding events are intermittent, i.e. a fluctuating number of many controlled 
perturbations of the hammer is necessary to induce a single sliding event of 
the cylinder. Each constant perturbation on the chute is associated with a 
time unit. The number of impacts of the hammer to obtain each sliding was 
recorded and from this data we performed the statistical analysis discussed in 
this paper. In all experiments the inclination was in the interval 12° to 18°, 
with 16° <9c< 32°, and the reduced angle {0^ — 0) /O^ was in the interval 0.19 
to 0.37 (see below for details). The sliding experiments were polarized in the 
sense that the cylinders were always placed on the chute in a same axial ori- 
entation. Two types of experiments were made with the azimuthal orientation 
of the cylinder fixed or not. In the first case, rotations of the cylinder around 
its symmetry axis are blocked, and as a consequence the interacting surfaces 
in contact were the same during all the experiment. The sliding spatial reso- 
lution is 0.5mm, and the angular resolution for the angle of inclination is 0.1°. 
Only sliding events with a length equal or larger than a threshold of 1mm 
are recorded for statistical purposes. The average sliding length in the exper- 
iments varies typically in the interval of 5 to 10mm. The experiments were 
reahzed in a small sealed room at a temperature of 24 (±1)°C, and 59 (±3)% 
of humidity. Chute and cylinder operated free of lubrication and both surfaces 
were submitted to a periodic cleaning with ethanol and cotton wool. The ex- 
perimental data obtained can be classified in four different groups, as follows. 

Group Gl consists of 9 time series of sliding events obtained in 1996 for massive 
cylinders with fixed diameter (f) = 9.6mm, and 9 different lengths L(mm) = 
5, 10, 20, 50, 100, 200, 300, 500, and 1000. Here and in all the other groups 
mentioned below, a series of sliding events is planned to generate a total of 
100 shding events. Both chute and cyhnder have nominal deviations of the 
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Fig. 1. Schematic diagram of the experimental apparatus: the inchne (A) has a 
Lo = 1,500mm long groove (the cross section is shown in B) and is rigidly main- 
tained at an inclination 6 with the horizontal. The cylinders (C) have lengths varying 
from L = 5mm to L = 1000mm. The hammer of m = 175g or 75g (D) impacts the 
base of the incline with a velocity close to 50cms~^. See text for detail. 

profile from a smooth surface in the range between about O.lyum and l/im 
[16]. The mass of the hammer used to introduce the perturbations on the 
rigid V-shaped chute was mi = 175g. The total duration of these time series 
varies from 364 (for L — 1000mm) to 6, 979 (for L — 5mm) time units or 
hammer impacts. In this group the azimuthal orientation of the cylinder is 
held fixed. The average (on the nine values of L previously mentioned) of the 
reduced angle of inclination in this group is {9^) = {{9c — 9) /9c) = 0.26 ±0.07. 

Group G2 consists of 27 time series obtained with the same set of cylinders 
and chute of Gl (here there are 3 independent time series associated with each 
one of the 9 values of L), but the data were collected in 2002, i.e. this group 
incorporates complex aging effects as for example the oxidation of the several 
interacting surfaces. Two relevant aspects in this group are: (i) before the 
beginning of the experiments, the aluminum cylinders and the chute had to 
be cleaned firstly with a fine steel wool, and after with water and neutral soap, 
(ii) The mass of the hammer used in the experiments in G2 was m2 — 75g. 
It is expected on the basis of these two modifications that the signal to noise 
ratio is somewhat inferior for G2 as compared with Gl. The effect of these 
modifications will be discussed in Sec. 3. Here, differently from Gl, both chute 
and cylinder have nominal deviations of the profile from a smooth surface in 
the range between about l//m and 10/xm [16] after the use of the steel wool. 
In this group, the experiments maintained the same interface cylinder-chute 
by preventing azimuthal rotations of the cylinder around its longitudinal axis. 
The average of the reduced angle of inclination in this group is the same as in 
Gl: {9,) = 0.26 ±0.07. 
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Group G3 refers to 17 time series of sliding events for the same cylinders 
and chute of G2, but the length L of the cylinders were restricted to the 
values L(mm) = 100, 200, 300, 500, and 1000. In the experiments grouped 
in G3, the control to maintain the cylinders free of rotations around their 
longitudinal axes was relaxed, i.e. in this case the cylinder-chute interface 
can vary continuously. Especially the cylinders with the smaller lengths are 
more sensitive to execute small rotations around their longitudinal axes. The 
hammer used in G3 was the same used in G2. In this group Oy- (not an average) 
is fixed at 0.37. 

Group G4 consists of 9 time scries of sliding events obtained with a new set of 
aluminum cylinders and rigid aluminum chute (also 1500mm long), with the 
same type of surface finishing used to obtain the data in Gl. In this group 
we have used hollow cylinders of 10mm external diameter and with 1mm wall 
width. The cylinders are free to rotate around their symmetry axes as in G3, 
and the hammer is the same used to obtain the data in G2 and G3. The length 
of the cylinders used in this case and the associated number of equivalent time 
series were, respectively lOOOmm/three scries, 500mm/two series, 200mm/two 
series, and lOOmm/two series. In this group 9y- is fixed at 0.25 for all the values 
of L. 

Before to conclude this Section, we would like to justify briefly the choice of 
some parameters in our experiments. Six years after the acquisition of the 
data grouped in Gl, we decided to test the robustness of this type of sliding 
experiment. As a consequence, the experiment with G2 was planned to include 
more time scries: performing three experiments for each value of the length L 
was considered a compromise between reliability and limitation of time. After 
a long period to obtain the data grouped in G2 we were forced to restrict the 
variabihty of the length of the cylinders for G3 and G4 (from 9 values of L in 
Gl and G2 to 5 values of L in G3, and 4 values in G4) while maintaining the 
approximate number of time series (2-3) for each value of length. This decision 
was taken to avoid a pronounced abrasion of the chute. To compensate this 
limitation, we increased our control on the reduced angle of inclination. 



3 Results and discussion 



We have measured in the experiments the number of impacts of the hammer, 
Tjv(I/), necessary to obtain N sliding events with size A > 1mm, for a cylinder 
with length L. To quantify the propensity or susceptibility of the cylinder 
to slide on the chute we introduce the sliding susceptibility Xn{L) deflned 
as the ratio of the number of sliding events or failures to stick, N, to the 
corresponding number of trials to move the cylinder through the controlled 
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perturbations on the chute, T]sr{L), that is 
Xn{L) = N/Tn{L). 



(1) 



Thus, Eq. (1) is an average quantity given by the reciprocal of the number 
of perturbations necessary to obtain one shding event. If a large number of 
perturbations is necessary to obtain only a small number of sliding events, the 
corresponding susceptibility will be small; on the contrary, if each perturbation 
of the hammer is followed by a sliding event, the susceptibility reaches its 
maximum value Xmax = 1- 

Figure 2 shows Xn{L) for = 50, 70, 85, and 100, for nine values of L in the 
interval 5 < L(mm) < 1000, for the data coming from the group Gl described 
in Section 2. Xn{L) in this case increases as the power-law L"^i(^), along 2.3 
decades of variability in L. The exponent cti ( A^) as a function of N, obtained 
after a best fit analysis of the experimental data for N varying from 1 to 100 
is given in the inset of Fig. 2. Its asymptotic value seems to be q;i(cxd) = 0.56. 
From Fig. 2 we can conclude that the scaling relation for Xn{L) is reasonably 
robust and independent of within typical statistical fluctuations of 5% to 
10% in the exponent a. Our overall estimate for this exponent from the data 
in Fig. 2 is q; = 0.55 ± 0.05. For N — 100, the average number of perturba- 
tions necessary for a single sliding event of the 5mm cylinder is close to 70 
(the experimental point associated with L = 5mm in Fig. 2 corresponds to 
an average on just 6,979 perturbations), whereas for the longest cylinder of 
1000mm, the corresponding average number of perturbations is close to 3.6 
(the experimental point associated with L — 1000mm in Fig. 2 corresponds to 
an average on just 364 perturbations). From now on we will focus on Xn{L) 
for N = 100, because for this particular value of the susceptibility function 
appears to be sufficiently close to its thermodynamic limit. Indeed, a detailed 
analysis (not shown in Fig. 2) indicates that the values of Xn{L) for L fixed 
seem to approach a hmit as A^ increases from 50 to 100. Moreover, for A^ = 100 
the exponent a also seems to have already reached its thermodynamic-limit 
value, as we can conclude from the inset of Fig. 2. For brevity we will then 
use the notation Xn=ioo{L) = xiL)- 

In Fig. 3 it is shown x{L) for G2 together with a power-law fit (continuous line) 
to the experimental data; here x{L) ~ L"'', with a2 = 0.33 ±0.04. Each point 
in Fig. 3 for a fixed value of L represents an average on three independent time 
series. In this case, the reduction of near 60% in the mass of the hammer as 
compared with Gl is accompanied by a reduction close to 40% in the scahng 
exponent. This reduction in the impact energy, as well as the alterations in 
the surfaces introduced by the cleanness of the cylinders and chute with steel 
wool and the aging effects introduced along the interval of six years between 
the experiments grouped in Gl and G2 (See Sec. 2) lead to a reduction in the 
signal/noise relation which is apparently materialized in the reduction of the 
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Fig. 2. Log-log plot of the experimental sliding susceptibility xn{L), for N = 50, 
70, 85, and 100, and 5 < L(mm) < 1000, for group Gl defined in Section 2. The 
straight line in the main figure represent the power-law best fit associated with 
these particular values of N. Our overall estimate from these data is the scaling 
relation x ~ 2^o.55±o.05_ inset shows the exponent ai(A^) obtained from the 
adjust XN ^ L'^'^W ^ for N varying from 1 to 100. The asymptotic value of the 
scaling exponent seems to be ai(cxD) = 0.56. See text, Sec. 3 for detail. 



exponent a and in the fluctuations in the power-law scaling (in particular for 
L = 100mm). 

The experimental dependence of x{L) for G3 is shown in Fig. 4; this suscep- 
tibility increases as L"^, with = 0.67 ± 0.06 for 100 < L(mm)< 1000. Each 
point in Fig. 4 for a fixed value of L represents an average on two to four 
independent time series. This average is important to obtain a good scaling 
behavior through the elimination of fluctuations associated with the tendency 
of the shortest cylinders to rotate around their symmetry axes. 

For hollow cylinders (G4) free to rotate around the longitudinal axis as with 
G3 we have obtained the data presented in Fig. 5, with x{L) ~ L°"^, = 
0.48 ± 0.02. In this case, to find a power-law fit similar to those appearing in 
Figs. 2 to 4, we need again to perform an average on similar time series to 
attenuate the fluctuations due to rotations. As noticed in Sec. 2, a total of 
9 time series were used to get Fig. 5: average on 3 series for cylinders with 
L = 1000mm; average on 2 time series for L = 500mm; average on 2 series for 
L = 200mm, and average on 2 series, for L = 100mm. 
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L (mm) 

Fig. 3. The same as in Fig. 2 but for group G2: the power-law fit to the experimental 
data gives x{^) ^ ^0-33±o.04_ g^^]^ point in this figure represents an average on three 
independent time series. See text, Sees. 2 and 3 for detail. 




L (mm) 

Fig. 4. Plot of the sliding susceptibility for G3; the best fit to the data gives 
X ~ ^o.67±o.06^ £qj. j^qq ^ L(mm) < 1000. As in Fig. 3, each point represents an 
average on independent time series. See text, Sees. 2 and 3 for detail. 
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L (mm) 

Fig. 5. Plot of the sliding susceptibility for hollow cylinders free to rotate around 

the longitudinal axis (group G4). In this case x{L) ~ L^-'^^^o.02_ ^igs. 3 and 4, 

each point in this plot represents an average on independent time series. See text, 
Sees. 2 and 3 for detail. 

The first physical aspect that we consider to explain the experimental results 
reported in this work is the effect of the vibrations induced on the chute and on 
the cylinder as a consequence of the impacts of the hammer. A simple heuristic 
argument predicts that the susceptibility Eq. (1) which is proportional to the 
probability of the cylinder to slip after a perturbation would scale linearly with 
the total number u of longitudinal/transversal/torsional vibration modes of 
the cylinder: 

X-^iy. (2) 

This is a reasonable assumption because (i) cylinders with more vibration 
modes tend to disengage more easily from the contacts that tend to maintain 
them fixed with respect to the chute, and (ii) the role of the chute on Eq. 
(2) can be eliminated because it is the same for cy finders of different lengths. 
On the other hand, for a thin rod-like system uniform along its length L we 
expect that 

i^-(VAmin)"°~(V0)"", a^ao = l, (3) 

where Amin — (= diameter of the rod/cylinder) is the minimum (cutoff) 
wavelength allowed [17]. If Eqs. (2) and (3) are assumed to be true, we should 
expect ;\; ~ 1/ ~ L"°, for constant; however, our data indicate a significantly 
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reduced value for the exponent a: the experimental values of a in the four 
groups of data examined in this work are of the order of 0.5 to 0.7aQ, 
for Gl and G3, to O.Sao to O.Sao, for G2 and G4. Our overall estimate is 
that the exponent a in this type of experiment can assume values in the 
interval 1/3 to 2/3, irrespective the group of data considered. We note in 
passing that Eqs. (2) and (3) imply that x would be independent of L if 
(j) ^ L, i.e. if the aspect ratio of the cylinder is held fixed, a result which can 
be experimentally tested. An important observation on Eq. (3) is that it is 
expected to be true only for simple boundary conditions, for instance, for a 
rod with fixed extremities or a rod with free extremities on a perfectly smooth 
non-interacting plane. It is natural to consider how the number of modes are 
modified by the roughness that is present in most real surfaces. For a rough 
cylinder of length L interacting with a rough fractal plane, v can scale as 
L°, where a is a non-trivial exponent smaller than the unit. This can occur, 
because in this case a fraction of the vibration modes will be forbidden or 
locked as a consequence of geometric engagements associated with the real 
topography of the interface. 

Alternately, Figs. 2 to 5 indicate that the average time between sliding events 
for a cyhnder of length L, r = Tn[L)/N — decreases as the power law 
T ~ L~°, with a varying typically from 1/3 to 2/3. It is interesting to notice 
that this behavior is reminiscent from the observed dependence of the average 
interval between earthquakes on a fault with length L: ~ L~^'^ [18]. Thus 
the interval between earthquakes on a specified fault is longer for smaller faults, 
and decays as a power-law, as similarly observed in our experiments for sliding 
events of rough solid cylinders on rough inclines. 

Now we will introduce in the discussion what seems to be the second physical 
aspect to be considered to explain the experimental results reported in this 
work: the effect of the roughness of the cyhnder-chute interacting surfaces. 
We conjecture that the origin of the tendency of x(^) for interacting rough 
surfaces to display power laws in L in several physical situations as exemplified 
in Figs. 2 to 5 is connected with the frequently observed fractal nonstationary 
behavior of the roughness of surfaces [1,15,19]; by this we understand that a 
sample of finite length L taken from a real surface will never, however long, 
completely represents its properties. If the height h{x) of a surface as a function 
of the position along a particular cut, x, is measured, the associated roughness 

1 /2 

can be defined by the width w = (^(^h{x)'^'j — {h(x))^^ , where {h(x)) — 

{1/ L) Jq h{x) dx, and (h{x)'^^ — h(xy dx. A great deal of experimental 
data [19,20,21,22] corroborate a robust scaling relation for w{L) along an 
impressive interval of variability of eight decades in L, from typically 10~®m 
to lO^m: 

w ~ ~ L^-^, (4) 
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where ( — 3 — D is the so caUed roughness exponent, and D is the surface 
fractal dimension in physical space [15,21]. If we use the scaling hypothesis 

xiL)r^w{L), (5) 



which says that the facility to generate sliding events as measured by x is pro- 
portional to the degree of roughness of the interacting fractal nonstationary 
surfaces (low(high) roughness =^ low(high) sliding susceptibility), the expo- 
nent a obtained in our experiments must be identified with the roughness 
exponent (. The reader can notice that Eq. (5) does not work for Euclidean 
rough surfaces, as for instance one with a regular saw tooth profile. In this 
last case, the roughness is stationary, w is independent of L (i.e. ( = 0), and 
as w grows, x is expected to decrease. Experimental data indicates that real 
surfaces have C ~ 1/2 for a multitude of natural and artificial surfaces, includ- 
ing concrete runways, hip joints, ship hull plates, and many types of machined 
metal surfaces [1,15,19,20], a result compatible with the exponent a obtained 
for Gl and G4 (see Figs. 2 and 5). Furthermore, the roughness of ductile cracks 
of aluminum alloys that received different heat treatments is also described 
by Eq. (4) with a roughness exponent about ( — 0.8 [21]. Another experi- 
mental study of the roughness of brittle cracks of Al-Si alloys, steel cooled 
by liquid nitrogen, graphite, porcelain, and other materials equally indicates 
the validity of Eq. (4) with ( close to 0.8 [22]. From Eqs. (4) and (5), and 
from the numerical values of the susceptibility exponent a(<-> Q obtained in 
Figs. 2 to 5, we conclude that the fractal dimension of the chute-cylinder in- 
teracting surfaces in our experiments are possibly random Brownian surfaces, 
whose fractal dimension is = 2.5, in agreement with the data reported in 
Ref. [19]. Why indeed should the sliding susceptibility be proportional to the 
roughness of the surface? There are in principle three arguments in favor of 
the conjecture given in Eq. (5): 

First of all, the susceptibility is according to Eq. (1) an average on fiuctuating 
sliding events. These fluctuations reflect obviously the fluctuating forces on 
the cylinders, which are in turn contact forces operating in the fractal nonsta- 
tionary interface cyhnder-chute. We should expect that these contact forces 
are affected by basic geometric aspects of the surfaces, and in particular by 
the surface roughness. Secondly, Eq. (5) reduces to the correct limit x = 0, 
for a perfect surface {w = 0). This is expected because in this case the chute- 
cylinder interacting surfaces are firmly welded in a single structure due to 
strong atomic interactions. Thirdly, if we consider what is said in the previous 
paragraph, and apply Eqs. (4) and (5) for earthquakes, we obtain x ^ L^'^ if 
we use the reported value D = 2.6 for the fractal dimension of fault surfaces 
[18]. The last scaling is equivalent to say that the predicted average time be- 
tween earthquakes on a fault with length L would scale as ~ L~° '', which 
is close to the reported value Te ~ L~°-^, obtained from seismic data [18]. 
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Moreover, if we use consistently Eqs. (2). (4) and (5), we obtain that the 
number of vibration modes accessible to the cylinder on the chute scales as 

V ~ L=^-^, (6) 

a result that reduces indeed to Eq. (3) only in the Euclidean limit D — > 2. In 
principle, we expect in conformity with Eq. (6), that cylinder and chute rough 
surfaces with increasing values of D tend to be more tightly engaged as a con- 
sequence of the increment in the effective area for microscopic interactions. 
This effect materializes itself in a microscopic reduction in the number of al- 
lowed vibration modes of the cylinder, i.e. v scales with L with an exponent 
necessarily smaller than the unit. In particular, a reduction in the number of 
allowed vibration modes is expected to be maximum if the interacting rough 
surfaces are space-filling STirfaces with sm^face fractal dimension D ^ 3. In 
this case, the surface interactions are saturated in the sense that they behave 
as volume-interactions that constrain the number of vibration modes to as- 
sume small values that are essentially independent of L (or per harps have a 
logarithmic dependence on the length of the cylinder) . For Brownian surfaces 
[p = 5/2), we obtain the result ^ ~ ~ L-^/^, which seems to describe the 
body of our experimental data. 



4 Conclusions 

We have introduced a statistical response function x that measures the de- 
gree of propensity of a rough cylinder to slide on a perturbed rough surface 
inclined below the angle of repose. This function was investigated from the 
experimental point of view and we have found that x increases with the length 
of the cylinder as a nontrivial power-law in a number of situations. A possible 
origin for this scaling behavior is discussed in terms of a connection with the 
number of vibration modes of a rough cylinder on a rough chute, as well as in 
terms of the statistics of the roughness of the interacting surfaces. 

It is natural to ask on the influence of the several parameters appearing in the 
experiments on the exponent a. Although a final answer to this subject can 
not be advanced at the present stage, some considerations can be outlined: (i) 
In the light of Eqs. (4) and (5), the reduction observed from cti = 0.55 ± 0.05 
to a2 — 0.33 ± 0.04 could be attributed to an effective increase in the surface 
fractal dimension of the cyhnders after the use of the steel wool as described 
in Sec. 2. (ii) The increase observed from a2 = 0.33 ±0.04 to = 0.67 ±0.06 
could be associated to the corresponding increase in the reduced angle of incli- 
nation from 9r = 0.26 to 9^ = 0.37. It can be noticed that the largest exponent 
(eta) is obtained in the experiment with the largest 9,.. (iii) The proximity be- 
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tween cci (massive cylinder) and 0:4 (hollow cylinder), both with almost the 
same 6,., could suggest that the response of hollow cylinders perturbed by a 
light hammer is equivalent to that of massive cylinders perturbed by a heavy 
hammer. However, due to the very time consuming nature of the experiments 
described in this paper, a number of experimental aspects remain to be sys- 
tematically investigated in the future. Thus, new experiments exploring the 
variation of the intensity of the perturbation on the chute and the variation 
of the materials involved, as well as the effect of the diameter of the cylin- 
ders are necessary to clarify this fluctuation phenomenon associated with the 
stick-slip dynamics. Furthermore, other important aspects to be answered are 
the type of variation of x with the reduced angle (^c — d)/&c, for L fixed; and 
the possibility of a universal value a = 1/2 for the susceptibility exponent in 
any experiment involving sliding of rough solids of the type discussed in this 
paper. 
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